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O ■ Abstract 
O 
(N 

£>Y A 'forward walking' Green's Function Monte Carlo algorithm is used to obtain 

expectation values for SU(3) lattice Yang- Mills theory in (3+1) dimensions. 

The ground state energy and Wilson loops are calculated, and the finite- 

size scaling behaviour is explored. Crude estimates of the string tension are 

derived, which agree with previous results at intermediate couplings; but more 

accurate results for larger loops will be required to establish scaling behaviour 

' at weak coupling. 

O ' 
O 

*Q '. I. INTRODUCTION 

o 
o 

Classical Monte Carlo simulations provide a very powerful and accurate method for the 
study of Euclidean lattice gauge theories. In the Hamiltonian formulation [[|], on the other 
hand, the corresponding quantum Monte Carlo methods have been somewhat neglected. 
^h. Here we present a study of SU(3) Yang-Mills theory in (3+1) dimensions, using the Green's 
Function Monte Carlo approach 0, adapted to lattice gauge theory by Chin et al ||. 

Quantum Monte Carlo methods in Hamiltonian lattice gauge theory have a somewhat 
chequered history. The first calculations used a strong-coupling basis involving discrete 



c3 



X 

"electric field" link variables, and a "Projector Monte Carlo" approach [HH, which used the 
Hamiltonian itself to project out the ground state. A later version of this was the "stochastic 
truncation" approach of Allton et al ||. Using this approach one can successfully compute 
string tensions and mass gaps for Abelian models |7]]. For non-Abelian models, however, 
some technical problems arose . The use of a Robson- Webber recoupling scheme |§ at the 
lattice vertices requires the use of Clebsch-Gordan coefficients or 6j-symbols, which are not 
known to high order for SU(3); and furthermore, the 'minus sign' problem rears its head, 
in that destructive interference occurs between different paths to the same final state. It 
may well be that a better choice of strong-coupling basis, such as the 'loop representation', 
might avoid these problems; but this has not yet been demonstrated. 

In the meantime, Heys and Stump and Chin et al. |J pioneered the use of "Green's 
Function Monte Carlo" (GFMC) or "Diffusion Monte Carlo" techniques in Hamiltonian 
LGT, in conjunction with a weak-coupling representation involving continuous gauge field 



link variables. This was successfully adapted to non-Abelian Yang-Mills theories |10Hl3 
with no minus sign problem arising. In this representation, however, one is simulating the 
wave function in gauge field configuration space by a discrete ensemble or density of random 
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walkers: it is not possible to determine the derivatives of the gauge fields for each config- 
uration, or to enforce Gauss's law explicitly, and the ensemble always relaxes back to the 
ground state sector. Hence one cannot compute the string tensions and mass gaps directly 
as Hamiltonian eigenvalues corresponding to ground states in different sectors, as one does 
in the strong-coupling representation. Instead, one is forced back to the more laborious 
approach used in Euclidean calculations: namely, to measure an appropriate correlation 
function, and estimate the mass gap as the inverse of the correlation length. We have in- 
troduced the 'forward- walking' technique, well-known in many-body theory PJT^-[T6|, to 



measure the expectation values and correlation functions. The technique has been demon- 
strated for the cases of the transverse Ising model in (1+1)D ||17| , and the U(l) LGT in 
(2+l)D 



Here we apply the technique for the first time to a non-Abelian model, namely SU(3) 
Yang-Mills theory in (3+l)D. The ground state energy and Wilson loop values are calculated, 
and approximate values are extracted for the string tension in the weak-coupling regime. 
Comparisons are made with earlier calculations, where they are available |19| . 



Our conclusions are that the method is a viable one, but requires the use of an improved 
"guiding wave function" to achieve better accuracy. There are certain drawbacks intrinsic 
to the method P,|2"0|], such as the necessity to use a branching algorithm and a guiding wave 
function, which tend to introduce substantial errors into the results, both statistical and 
systematic. For these reasons, it may be preferable to employ a Path Integral Monte Carlo 
approach to these models, which avoids the problems above. 

Our methods are presented in Section II, the results are outlined in Section III, and our 
conclusions are discussed in Section IV. 



II. METHOD 
A. Lattice Hamiltonian 



The Green's Function Monte Carlo formalism has been adapted to SU(2) Yang-Mills 
theory by Chin, van Roosmalen, Umland and Koonin JTT|, and sketched for the SU(3) case 



by Chin, Long and Robson [O. Here we provide a slightly fuller discussion of the SU(3) 



case, following the earlier treatment of Chin, van Roosmalen et al JTT 
The SU(3) lattice Hamiltonian is given by 

H = f{EW " \H T <U P + Ul)} (1) 
za I 6 p 

where is a component of the electric field at link 1, A = Q/g A , the index a runs over the 8 
generators of SU(3), and U p denotes the product of four link operators around an elementary 
plaquette. 

Commutation relations between electric field and link operators are 

[EtM^X-UA,, (2) 

choosing the Ef as left generators of SU(3), where the {A a ,a = 1,8} are the Gell-Mann 
matrices for SU(3). We will work with the dimensionless operator 
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H=W E i E i-^E Tr (U P + Ul) (3) 
z i p 

The link variables are elements of the group SU(3) in the fundamental representation 

U = exp(-i~AM a ) (4) 

There is no simple equivalent of the quaternion representation for SU(2). Following Beg and 
Ruegg pll ) we can represent 

' Zi z 2 z 3 \ 

U = ui u 2 u 3 (5) 
\ Wi w 2 w 3 J 

where z,u,w are three-dimensional complex vectors; then if U is to be unitary, we require 
z,u,w to be orthonormal; and if U is to have determinant unity, we require 

eijkZiUjWk = 1 (6) 

which is satisfied if 

Wi = e ijk z*u* k (7) 
One possible representation which satisfies these conditions is 

21 = (1 - i(x 3 + X & /y/3))/N 1} 

z 2 = {-x 2 - ix^/Ni, 

z 3 = (-x s —ix4 L )/N 1 (8) 



where 



and 



where 



and 



Ni = [l + (x 3 + x 8 /v/3) 2 + xj + x\ + x\ + x 2 5 } 1/2 ; (9) 

u 3 = (-x 7 - ix 6 )/N 2 = u 3 /N 2 , (10) 

u 2 = (1 - *(x 8/ V3 - ar 3 ))/JVi = u 2 /N 2 , (11) 

Mi = u~x/n 2 , ui = -\u 2 z* 2 + u 3 z 3 ]/zl (12) 

iv 2 = [N 2 + |« 2 | 2 + N 2 ] 1/2 (13) 



Wi = e ijk z*u* k . (14) 
This involves 8 unrestricted parameters {x a ,a = 1..8}, as expected. For small x, 

U^I-i\ a x a (15) 

i.e. 

x a w A a /2, (16) 

where A a is the gauge field on that link. The product of two link variables can be found by 
simple matrix multiplication. 
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B. Green's Function Monte Carlo method 



The Green's Function Monte Carlo method employs the operator exp(— t(H — E)), i.e. 
the time evolution operator in imaginary time, as a projector onto the ground state |"0o) : 

|Vo> oc \im{e- T{H - E) \m 

lim e -NAr(H-E)^ (17) 

where |$) is any suitable trial state. To procure some variational guidance, one performs a 
"similarity transformation" with the trial wave function $, and evolves the product ^\i/jq) 
in imaginary time. The heart of the procedure is the calculation of the matrix element 
corrersponding to a single small time step At. Chin et al [[l ij show that 



(x.'\<l>e- ATiH - E) <l>- 1 \x.) = Y[(U[\N{exp{--ATE?E?) exp[ArE?(Ef In ^)}}\U l ) 

i ^ 

exp{Ar[£ - ^iJ^x)]} + 0(Ar 2 ) 

= p(x',x)w(x) + 0(Ar 2 ) (18) 

where x = {U{\ denotes an entire lattice configuration of link fields. 

The product §\ip) is simulated by the density of an ensemble of random walkers, as in the 
SU(2) case. At the kth. step, the 'weight' of each walker at x^ is multiplied by u'(xfc) and 
the next ensemble {x^+i} is evolved from {x^} according to the matrix element p(xk+i, Xk). 
The effect of p(xk+i, Xk) is to alter each link variable Ui in {xk} to U( by a Gaussian random 
walk plus a "drift step" guided by the trial wave function: 

U' = AUU d U (19) 

where Ua = exp[i^X a (iArE a In $)] is the drift step, and AU is an SU(3) group element 
randomly chosen from a Gaussian distribution around the identity, with variance (As 2 ) = 
8 At (i.e. At for each index a), where 

(As 2 ) w ]T A a A a = 8Ar, (small A a ) (20) 

a 

or 

(x a x a ) = ^(A a A a ) « each a. (21) 

The simulation is carried out for a large number of iterations At, until an equilibrium 
distribution <^> |^ ) is reached. The energy E in (|ll|) is adjusted after each iteration so as to 
maintain the total ensemble weight constant. The average value of E can then be taken as 
an estimate of E , the ground-state energy. 

As time evolves, the weights of some walkers grow larger, while others grow smaller, 
which would produce an increased statistical error. To avoid this, a "branching" process is 
employed, whereby a walker with weight larger than some threshold is split into two inde- 
pendent walkers, while others with weights lower than another threshold are amalgamated. 
We use Runge's technique |22| for this purpose. 
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C. Trial Wave Function 



The trial wave function is chosen to be the one-parameter form [13 

$ = exp[aJ2Tr(U p + Ul)} 



(22) 



Then the drift step is |TT 



for each link, where 



U d = exp[iy(iAr£ a ln$)] 
= exp[-i^Af] 



(23) 



Af = -iAT-Y.TrlXV^Ut - h.c] 

1 pel 



(24) 



i.e. 



■T, 



Af 
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pel 



(25) 



(note: the effect of Ef on a plaquette operator is to 'insert' a A a /2 in front of the appropriate 
link operator C/j, to be followed by the remaining link operators in the plaquette, taken in 
the direction of the link 1). 

Finally, the trial energy factor is 

= y (£Tr[A a ^..[/| - h.c.]) 2 + (^ - A) Y,Tr{U p + U})}. (26) 
i pel pel 

Therefore the weight factor is 

w(x) = exp{Ar(£ - ^H®)} 

exp{Ar(E trml - _ A) £ £ 2Re{TrU p } 



24' 



i pel 



+ ^E(E2^{Tr[A^..^]]}) 2 ))}. 

i pel 



(27) 



D. Forward Walking estimates 

The "forward walking" technique is used to estimate expectation values M. Its appli- 



cation to the U(l) lattice gauge theory in (2+l)D was discussed by Hamer et al | 18| . It is 
based on the following equation: for an operator Q, 
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<Q>0 (MOM 
($ |# j q|^ > 



J — »oo 



}$|i^|Vo> 

E K(*j, xj_i) . . . K(x 2 , xi)Q(xi)^ (xi) 

E ^(Xj, Xj_i) . . . ^(x 2 , Xi)^o(Xi) 



(28) 



where i^(xj,xj_i) is the evolution operator for time At, and .fT(xj, xj_i) is the same 
operator in the similarity transformed basis. Again we have assumed that the operator Q 
is diagonal in the basis of plaquette variables x. 

This equation is implemented by [p^ -[T6[| the following procedure: 

i) Record the value Q(xj) for each "ancestor" walker at the beginning of a measurement; 

ii) Propagate the ensemble as normal for J iterations, keeping a record of the "ancestor" 
of each walker in the current population; 

hi) Take the weighted average of the Q(xj) with respect to the weights of the descendants 
of Xj after the J iterations, using sufficient iterations J that the estimate reaches a 
'plateau'. 



III. RESULTS 

Simulations were carried out for LxLxL lattices up to L=8 sites, using runs of typically 
4000 iterations and an ensemble size of 250 to 1000 depending (inversely) on lattice size. 
These statistics are approximately 100 times less than those used in the U(l) calculations 
|l~8|, but about 100 times greater than those used in the previous SU(3) calculations of 



Chin et al |13|]. Time steps At of 0.01 and 0.05 "seconds" were used, with each iteration 
consisting of 5 sweeps and 1 sweep through the lattice, respectivelyfollowed by a branching 
process. The first 400 iterations were discarded to allow for equilibration. The data were 
block averaged over blocks of up to 256 iterations, to minimize the effect of correlations on 
the error estimates. 

The results taken at At = 0.01 and At = 0.05 were extrapolated linearly to At = 0. 
Figure 1 demonstrates that the dependence of the ground-state energy on At is approxi- 
mately linear. 

The variational parameter c was given values as shown in Table I. These are essentially 



the values used by Chin et al ||13|| , obtained from a variational Monte Carlo calculation. We 
checked that these were approximately the optimum values for small lattices. 

Forward-walking measurements were taken over J iterations, where J ranged from 20 to 
100, depending on the coupling A. Ten separate measurements were taken over this time 
interval, in order to check whether the value measured by forward-walking had reached 
equilibrium. A new measurement was started soon after the previous one had finished. 
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A. Ground-state Energy 



The dependence of the ground-state energy on the variational parameter c is illustrated 
in Figure 2. It can be seen that the energy reaches a broad minimum at about the expected 
value (c = 0.33 at this coupling). 

Our estimates of the ground-state energy are listed in Table I, as a function of the 
coupling A and lattice size L. The dependence on lattice size is illustrated in Figure 3, at 
two fixed couplings A = 3.0 and A = 5.0. In the "strong-coupling" case, A = 3.0, it can 
be seen that the results converge exponentially fast in L, whereas in the "weak-coupling" 
regime, A = 5.0, the convergence is more like 1/L 4 at these lattice sizes. This behaviour 
merits some further explanation. 

A similar phenomenon occurs in the case of the U(l) theory in (2+l)D [^,^,|T^]. In 
the strong-coupling regime, where the mass gap is large, the usual exponential convergence 
occurs. In the weak-coupling regime, however, where the mass gap M is very small, the 
finite-size scaling behaviour for small lattice sizes is that of a massless theory, and it is only 
at much larger lattice sizes L ~ 1/M that a crossover to exponential convergence occurs. In 
the U(l) case, it has been shown [p|, |l8[| that the finite-size scaling behaviour at small L is 
well described by an "effective Lagrangian" approach, using the Lagrangian for free, massless 
photons that the model was originally constructed to simulate. In the same way, an "effective 
Lagrangian" corresponding to free, massless gluons (non-interacting QCD) should describe 
the finite-size behaviour in the present case, in line with the idea of asymptotic freedom. 
By analogy with the (2+l)D case, we expect a 1/L 4 dependence for the corrections to the 
ground-state energy per site. We hope to pursue this analysis further at a later date. 

An anomalous feature in Figure 3b) is that the L = 8 point lies well out of line with the 
others. This occurs at other couplings also. We suspect that the results for L = 8 are not 
reliable, and that the trial wave function will have to be further improved to give reliable 
results for such large lattices. Supporting evidence for this will be presented below. 

We have made estimates of the bulk limit, extrapolating mainly from the smaller L values 
where possible, and the results are listed in Table II. Our present estimates generally lie a 
little below those of Chin et al [Ol, and we believe them to be more accurate in view of 



our greater statistics. The estimates for the bulk ground-state energy per site are graphed 
as a function of coupling in Figure 4, where they are compared with previous estimates 
|i~9|| obtained by an 'Exact Linked Cluster Expansion' (ELCE) procedure, and with the 
asymptotic weak-coupling series pl| 



e 



-3A + 7.798A 1/2 , A -> oo. (29) 



The Monte Carlo results agree very well with the ELCE estimates, and appear to match 
nicely onto the expected weak-coupling behaviour for A > 6. 



B. Wilson Loops 

The forward-walking method was used to estimate values for the mxn Wilson loops, 
W(m,n). Figure 5 shows an example, namely the estimate of the mean plaquette W(l, 1) 
as a function of J for the case L = 6, A = 1.5. It can be seen that the estimate relaxes 
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exponentially towards a plateau value as the number of iterations J increases: an exponential 
fit is performed to estimate the asymptotic value. It can also be seen that the statistical 
error on each point is much larger than the point-to-point variation: the estimates are 
highly correlated, and all the points tend to move up and down together as one goes from 
one sample to the next. 

A problem that arises in these measurements is the loss in statistical accuracy at large 
couplings. At large couplings the weights of the random walkers vary rapidly with time, and 
it can easily happen that during a measurement the descendants of all the 'ancestor' walkers 
but one die out from the ensemble, at which point the result 'freezes', and the number of 
members of the ensemble has effectively been reduced to one. This inevitably means a severe 
loss of statistical accuracy. 

A graph of the 'mean plaquette' W(l, 1) versus the variational parameter c is shown 
in Figure 6. Another problem is immediately apparent. The estimate for W(l, 1) is not 
independent of c, in fact it depends linearly on c over this range, and the size of the variation 
is such that the probable systematic error due to the choice of c is an order of magnitude 
larger than the random statistical error in the results. Thus it would be advatageous in 
future studies to put more effort into improving the trial wave function, rather than merely 
improving the statistics. 

Figure 7 shows examples of the dependence of the results on lattice size L. Once again, 
the results at strong coupling A = 1.5 converge exponentially fast, while those at the weak 
coupling value A = 5.0 can be approximately fitted by a 1/L 4 dependence. The L = 8 value 
and even the L = 6 value again lie off the trend of the smaller lattices, and are probably not 
very reliable. 

Estimates of the bulk limit are listed in Table II. The estimates for the mean plaquette 
are graphed as a function of coupling A in Figure 8, and compared with series estimates at 
strong and weak coupling |l9|,p4| . The agreement is quite good. 



C. String Tension 



Having obtained estimates for the Wilson loop values on the bulk lattice, one can extract 
estimates for the 'spacelike' string tension using the Creutz ratios: 



Ka 2 ~ R n = - In 



W(n, n)W(n — 1, n — 1) 
W(n,n- l) 2 



(30) 



or the cruder 2-point estimates 



-I In 

n 



W(n,n) 
W(n,n - 1) 



(31) 



The results are shown in Figure 9. Also shown in Figure 9 are some previous estimates 
derived from the 'axial' string tension, obtained |H5|] using an 'Exact Linked Cluster Expan- 
sion' (ELCE) method. The axial string tension aT is calculated as an energy per link, and 
must be converted to a dimensionless, 'spacelike' tension by dividing by the 'speed of light' 
c, 



S 



c 



(32) 



where 26 



c ~ 



0.1671s 2 




l-0.1671i 



oo. 



(33) 



We have also used the weak-coupling relationship between the scales of Euclidean and Hamil- 
tonian lattice Yang-Mills theory calculated by Hasenfratz et al [^] to plot the results against 
the Euclidean coupling (3 = Q/g\ 



f E) where 



6A- 0.01308 



(34) 



It can be seen that the present GFMC results are in rough agreement with the axial 
string tension results in the region 4 < (3 < 5, which is also the region where the 'roughening' 
transition occurs in the string tension [IS]. For (3 > 5, however, the Creutz ratio R2 runs 
above the ELCE estimate, and shows no sign of the expected crossover to an exponentially 
decreasing scaling behaviour at (3 ~ 6. We presume that this is a finite- size effect, and that 
the Creutz ratios R n for larger n will show a substantial decrease in the 'weak-coupling' 
regime f3 > 6. That is certainly the pattern seen in the Euclidean calculations ||25|| , or in the 



U(l)2+i model [18|. Unfortunately, however, our present results for the larger Wilson loops 
are not of sufficient accuracy to allow worthwhile estimates of R n for n > 2. 



IV. SUMMARY AND CONCLUSIONS 

We have presented the results of a new Green's Function Monte Carlo study of the SU(3) 
Yang-Mills theory in the (3+l)D Hamiltonian formulation. A forward- walking method has 
been used to estimate values for the Wilson loops as well as the ground-state energy, and 
hence some rather crude estimates of the string tension have been extracted. Comparisons 
have been made with an earlier Hamiltonian calculation of the axial string tension ||19| . The 
two sets of results agree in the 'roughening' region; but our Monte Carlo results do not 
extend to the large Wilson loops that would be required to demonstrate 'scaling' behaviour 
in the weak-coupling regime. 

Some significant problems with the GFMC method have emerged from this study. The 
'forward-walking' technique was introduced specifically to avoid any variational bias from 
the trial wave function PJl4]-4T6|- As it turns out, however, the results for the Wilson loops 
show a substantial dependence on the trial wave function parameter c. The systematic error 
due to this dependence is an order of magnitude larger than the statistical error, so it would 
pay to put more effort in future studies into improving the trial wave function, rather than 
simply increasing the statistics. Furthermore, the effective ensemble size decreases during 
each measurement as the descendants of each 'ancestor' state die out, and this produces a 
substantial loss in statistical accuracy at weak coupling, as well. 

It would be preferable if one were able to do away entirely with all the paraphernalia 
of trial wave function, weights, branching algorithms, etc, and just rely on some sort of 
Metropolis-style accept/reject algorithm to produce a correct distribution of walkers. Within 
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a quantum Hamiltonian framework, a way is known to do this, namely the Path Integral 
Monte Carlo (PIMC) approach [pOR . We conclude that the PIMC approach may be better 
suited than GFMC to the study of large and complicated lattice Hamiltonian systems. 
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TABLES 



A 


1.5 


3.0 


4.0 


5.0 


6.0 


9.0 


c 


0.10 


0.20 


0.27 


0.33 


0.40 


0.58 


Ground state energy per site, e 


L = 2 


-0.15449(2) 


-0.68081(1) 


-1.3212(5) 


-2.359(1) 


-3.666(2) 


-8.389(6) 


3 


-0.15438(4) 


-0/6764(2) 


-1.2846(3) 


-2.187(1) 


-3.431(2) 


-8.042(4) 


4 


-0.15436(1) 


-0.6764(1) 


-1.2836(8) 


-2.160(3) 


-3.322(4) 


-7.94(1) 


6 


-0.15436(2) 


-0.6765(1) 


-1.2809(7) 


-2.133(2) 


-3.234(3) 


-7.613(9) 


8 




-0.6759(2) 


-1.2705(4) 


-2.094(2) 


-3.154(2) 


-7.402(6) 


Wilson loops, W(l,l) 


L = 2 


0.0721(2) 


0.1671(3) 


0.273(2) 


0.396(6) 


0.460(7) 


0.576(4) 


3 


0.0718(2) 


0.1641(4) 


0.245(2) 


0.346(4) 


0.437(5) 


0.562(2) 


4 


0.0718(2) 


0.1647(7) 


0.240(2) 


0.330(4) 


0.401(5) 


0.556(4) 


6 


0.0718(1) 


0.1648(4) 


0.240(3) 


0.310(1) 


0.383(2) 


0.550(3) 


8 




0.1633(8) 


0.228(1) 


0.296(1) 


0.365(2) 


0.545(2) 


W(2,l) 


3 


0.0059(1) 


0.0318(3) 


0.070(1) 


0.146(4) 


0.234(6) 


0.376(3) 


4 


0.0061(1) 


0.0322(4) 


0.067(2) 


0.130(4) 


0.182(4) 


0.355(6) 


6 


0.0061(1) 


0.0322(3) 


0.069(2) 


0.109(1) 


0.165(2) 


0.329(3) 


8 




0.0309(5) 


0.057(1) 


0.097(1) 


0.144(2) 


0.313(2) 


W(2,2) 


4 




0.0010(7) 


0.006(2) 


0.022(5) 


0.044(5) 


0.183(7) 


6 




0.0022(3) 


0.009(2) 


0.018(1) 


0.037(2) 


0.137(4) 


8 




0.0011(6) 


0.005(1) 


0.013(1) 


0.023(2) 


0.111(3) 



W(3,2) 



6 








0.004(1) 


0.006(2) 


0.059(3) 


8 








0.001(1) 


0.004(1) 


0.039(2) 



W(3,3) 



6 












0.020(3) 


8 












0.009(2) 



TABLE I. Variational parameter c, ground state energy per site e a and Wilson loops W(n, m) 
as functions of the lattice size L and coupling A. 



A 


1.5 


3.0 


4.0 


5.0 


6.0 


9.0 


Ground state energy per site, e 


this work 
Chin et al [13] 


-0.15436(2) 


-0.6764(1) 
-0.675(0) 


-1.284(2) 
-1.275(3) 


-2.16(2) 
-2.088(6) 


-3.25(2) 
-3.183(6) 


-7.9(1) 
-7.50(3) 


Wilson loops W(l,l) 


this work 
Chin et al [13] 


0.0718(2) 


0.165(1) 
0.1605(5) 


0.240(1) 


0.32(1) 
0.298(5) 


0.39(1) 
0.377(3) 


0.554(3) 
0.539(3) 


W(2,l) 


this work 


0.0061(2) 


0.0322(5) 


0.068(2) 


0.12(1) 


0.165(5) 


0.34(1) 


W(2,2) 


this work 




0.0021(1) 


0.008(4) 


0.018(4) 


0.035(7) 


0.13(1) 


W(3,2) 


this work 








0.004(2) 


0.006(2) 


0.05(1) 


W(3,3) 


this work 












0.01(1) 



TABLE II. Estimates of the bulk ground-state energy per site and Wilson loops as functions 
of A. 
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FIGURES 



FIG. 1. Estimated ground-state energy for lattice size L = 4 and coupling A = 3.0, c = 0.20, 
as a function of time step At. 

FIG. 2. Estimated ground-state energy as a function of the variational parameter c, for 
L = 6, A = 5.0. 

FIG. 3. Ground-state energy per site graphed against 1/L 4 , where L is the lattice size: a) at 
coupling A = 3.0, b) at coupling A = 5.0. The lines are merely to guide the eye. 

FIG. 4. The bulk ground-state energy per site graphed against coupling A. The points are our 
Monte Carlo estimates; the solid line represents earlier ELCE estimates [19]; and the dashed line 
represents the asymptotic weak-coupling behaviour. 

FIG. 5. Measured value for the Wilson loop W(l, 1) as a function of the number of for- 
ward-walking iterations J, for the case L = 6, A = 1.5. The solid line is an exponential fit to the 
data. 

FIG. 6. Estimated value for the mean plaquette W(l, 1) as a function of the variational 
parameter c, for L = 6, A = 5.0. 

FIG. 7. The mean plaquette W(l, 1) graphed against 1/L 4 , where L is the lattice size: a) at 
coupling A = 1.5, b) at coupling A = 5.0. The lines are merely to guide the eye. 

FIG. 8. The mean plaquette W(l, 1) for the bulk system graphed against coupling A. The 
solid line represents the strong-coupling series expansion [19], and the dashed line the asymptotic 
weak-coupling behaviour. 

FIG. 9. The string tension Ka 2 graphed against coupling f3. The circles are obtained from 
ELCE estimates of the axial string tension[19]; the triangles are Monte Carlo estimates of i?2- 
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